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Abstract: We study a 4d supersymmetric matrix model with a cubic term, which incor- 
porates fuzzy spheres as classical solutions, using Monte Carlo simulations and perturbative 
calculations. The fuzzy sphere in the supersymmetric model turns out to be always stable 
if the large- iV limit is taken in such a way that various correlation functions scale. This 
is in striking contrast to analogous bosonic models, where the fuzzy sphere decays into 
the pure Yang- Mills vacuum due to quantum effects when the coefficient of the cubic term 
becomes smaller than a critical value. We also find that the power-law tail of the eigenvalue 
distribution, which exists in the supersymmetric model without the cubic term, disappears 
in the presence of the fuzzy sphere in the large- N limit. Coincident fuzzy spheres turn out 
to be unstable, which implies that the dynamically generated gauge group is U(l). 
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1. Introduction 

Matrix models are considered as one of the most promising candidates for a nonperturbative 
formulation of string theories. Indeed some concrete models are proposed as constructive 
definitions of superstring and M theories [1,2]. These models are obtained from the di- 
mensional reduction of super Yang-Mills theory in ten dimensions. In a broad sense, such 
models belong to the class of the so-called dimensionally reduced models [3] (or large-iV 
reduced models), which were studied intensively in the eighties as an equivalent description 
of large-iV gauge theories. Unlike the old models, however, the new models are written in 
terms of Hermitian matrices, and they have manifest supersymmetry, which is expected to 
have crucial effects on their dynamics 1 . 

1 See ref. [4] for a comprehensive review on these issues. 
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An important feature of these matrix models is that the space-time is not introduced 
from the outset, but it emerges dynamically as the eigenvalue distribution of the bosonic 
matrices. In fact there are certain evidences in the IIB matrix model [2] that four- 
dimensional space-time is generated dynamically [5-8]. In refs. [5-7] the free energy of 
space-time with various dimensionality has been calculated using the gaussian expansion 
method, and the free energy turned out to take the minimum value for the four-dimensional 
space-time. In ref. [8] it was found that the fuzzy S 2 x S 2 (but not the fuzzy S 2 ) is a solution 
to the 2-loop effective action. See also refs. [9-20] for related works on this issue. 

By adding a Chern-Simons term to the matrix models, one obtains fuzzy spheres [21] as 
classical solutions [22-26], and their dynamical properties have been studied in refs. [27-39]. 
This provides a matrix description of the so-called Myers effect in string theory [40] . The 
emergence of a fuzzy sphere in matrix models may be regarded as a prototype of the 
dynamical generation of space-time since it has lower dimensionality than the original 
dimensionality that the model can actually describe. When k fuzzy spheres coincide, the 
gauge symmetry enhances from U(l) fc to U(k). By expanding the theory around such a 
solution, one obtains a U(fc) gauge theory on a noncommutative geometry [23]. Therefore 
the model may also serve as a toy model for the dynamical generation of gauge group, 
which is expected to occur in the IIB matrix model [41]. 

In fact one can use the above matrix models to define a regularized field theory on the 
fuzzy sphere as has been done on a noncommutative torus [42], which enables nonpertur- 
bative studies of such theories from first principles [43] . This is motivated from the general 
expectation that noncommutative geometry provides a crucial link to string theory [44] 
and quantum gravity [45]. Yet another motivation is to use the fuzzy sphere (or its gen- 
eralization [46,47]) as a regularization scheme alternative to the lattice regularization [48]. 
Unlike the lattice, fuzzy spheres preserve the continuous symmetries of the space-time con- 
sidered, and hence it is expected to ameliorate the well-known problem concerning chiral 
symmetry [36, 49-63] and supersymmetry. A challenge in this direction is to remove the 
effects of noncommutativity of the space-time in the "continuum limit" . The fuzzy sphere 
is also useful in the Coset Space Dimensional Reduction [64,65], where one can take the 
compact part of space-time to be a fuzzy coset [66,67]. 

Whatever the motivation is, the stability of fuzzy-sphere-like solutions is clearly one 
of the most important issues. In cases when there are more than one stable solutions, one 
can identify the true vacuum by comparing the corresponding free energy. This will be 
important in the dynamical determination of the space-time dimensionality and the gauge 
group in superstring theory. In the series of papers [68-73], we addressed such issues in 
various kinds of models using both perturbative calculations and Monte Carlo simulations. 
In ref. [68] we have studied the dimensionally reduced 3d Yang-Mills models with the 
Chern-Simons term, which has the fuzzy 2-sphere (S 2 ) as a classical solution [23]. We have 
found a first-order phase transition as we vary the coefficient (a) of the Chern-Simons term. 
For small a the large- ./V behavior of the model is the same as in the pure Yang-Mills model, 
whereas for large a a single fuzzy S 2 appears dynamically. In addition we find that the k 
coincident fuzzy spheres, which are also classical solutions of the same model, cannot be 
realized as the true vacuum in this model even in the large- N limit. This implies that the 
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dynamical gauge group is U(l) in this model. In refs. [69,70,73] we have extended this work 
to various matrix models, which incorporate four-dimensional fuzzy manifolds as classical 
solutions. While the fuzzy S 4 turned out to be unstable [69], we find that the fuzzy CP 2 [70] 
and the fuzzy S 2 x S 2 [73] are stable at large N although the true vacuum is actually given 
by the fuzzy S 2 . In the latter two cases the gauge group generated dynamically turned out 
to be U(l) as well. In ref. [72], on the other hand, it has been shown for the first time that 
gauge groups of higher rank can be realized in the true vacuum by adding a mass term to 
the 3d Yang-Mills-Chern-Simons model. 

The aim of the present paper is to study the impact of supersymmetry on the dynamics 
of the fuzzy spheres. The simplest 3d supersymmetric model [23] is problematic nonper- 
turbatively since the partition function is divergent [74-78] . This leads us to study the 4d 
supersymmetric model with a cubic term instead. Indeed it turns out that the supersym- 
metry has a striking impact. Unlike the bosonic models, the fuzzy sphere is always stable 
if the large- N limit is taken in such a way that various correlation functions scale. We also 
observe an interesting phenomenon that the power-law tail of the eigenvalue distribution, 
which exists in the supersymmetric models without the Chern-Simons term [77,79], disap- 
pears in the presence of the fuzzy sphere in the \arge-N limit. Coincident fuzzy spheres 
turn out to be unstable, which implies that the dynamically generated gauge group is U(l) 
in the present model. 

This paper is organized as follows. In section 2 we define the model and discuss its 
fuzzy sphere solutions. In section 3 we study the phase diagram of the model. In section 4 
we study the geometrical structure of the dominant configurations. In section 5 we study 
coincident fuzzy spheres and discuss the dynamical gauge group. Section 6 is devoted to 
summary and discussions. In appendix A we explain the algorithm for our Monte Carlo 
simulations. In appendices B and C we provide the details of perturbative calculations. 

2. The model and the fuzzy sphere 

The model we study in this paper is defined by the action 

S = S h + S { , (2.1) 

S h = Nti I -\[A^ A u ] 2 + ha ei jk AiAjA k \, (2.2) 

Si = -Nti^ a (T^ a/3 [A^ V/?]) , (2-3) 

where A^ (fx = 1,2,3,4) are N x ./V traceless hermitian (bosonic) matrices, and ipa^^Pa 
(a = 1,2) are N x N traceless complex (fermionic) matrices. Here and henceforth we assume 
that repeated Greek indices are summed over all possible integers. The k = 1, 2, 3) 

is a totally anti-symmetric tensor with ei23 = 1. The 2x2 matrices are Weyl-projected 
gamma matrices in four dimensions, and they are given explicitly as 

(2 - 4) 
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The convergence of the integration over is a non-trivial issue since the integration 
region is non-compact. At a = the partition function is finite for arbitrary N, as first 
conjectured by ref. [74] and proved later by ref. [77], and this remains to be the case also 
for a / [78]. Moreover, since the fermion determinant of this model is positive semi- 
definite (See ref. [11] for a proof), the model can be studied by Monte Carlo simulations 
without confronting the so-called sign problem. The pure super Yang-Mills model (a = 0), 
which may be regarded as the 4d version of the IIB matrix model [2], has been studied 
intensively [11,14,15,74]. The sign problem does not occur even if one includes the cubic 
term, which is real 2 . 

For a = the model has manifest SO(4) symmetry and super symmetry. The cubic 
term in (2.2) obviously breaks the SO(4) symmetry down to SO(3). It also breaks supersym- 
metry, but the effects of breaking is "soft" since the power of is lower than the quartic 
term [31]. Therefore one may still anticipate to see peculiar effects of supersymmetry. We 
repeat that the 3d supersymmetric model, which has been studied perturbatively [23,31], is 
actually problematic nonperturbatively since the partition function is divergent [74,77,78]. 
Therefore, the present 4d model is the simplest model that can be studied in order to 
examine the impact of supersymmetry on the fuzzy sphere dynamics. 

Let us then consider the classical solutions of this model. For V = the equation of 
motion reads 

3 

[A v , [A v ,Ai]] + ia e ijk [Aj,A k ] = for i = 1,2,3 , 
j,k=i 

[A v ,[A v ,At]] =0 . (2.5) 

Apart from the solution given by commuting matrices, which exists also for a = 0, we have 
the fuzzy S 2 solution given by 

Uf>=aL\ N) fori = 1,2,3, 

U sl, = o, <6) 

where l}p (i = 1, 2, 3) represents the r-dimensional irreducible representation of the SU(2) 
Lie algebra 

[4\L^] = iae ijk 4K (2.7) 

The solution A^ satisfies 

y t(Af ) f = \(N*-l)c?l If , (2.8) 
i=i 

which implies that the "radius" of the fuzzy sphere is given by 

p = -a VN 2 - 1 . (2.9) 

We consider more general solutions in section 5. 

2 This is in contrast to the Chern-Simon term in ordinary gauge theories in euclidean space-time, which 
is purely imaginary. Note that the coefficient a in (2.2) should be chosen to be real in order for fuzzy sphere 



solutions to exist. 
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3. Phase diagram 



3.1 Monte Carlo results 

In this section we calculate various quantities by Monte Carlo simulation and study the 
phase diagram of the model (2.1). We show results obtained by using the fuzzy sphere Ajj, 
as the initial configuration, but we have checked that the result is the same for other initial 
configurations such as = or some randomly generated configurations. For brevity we 
introduce the notation 



F\xv — i [Afi, A v \ , 
2 3 

M = Jn 1 Yl e ljk tv{AiA 3 A k ' 

i,j,k=l 



(3.1) 
(3.2) 



We note that there is an exact result 

1 



v tr(iVH+3a(M) = 6(l y2 



(3.3) 



which can be derived as in the bosonic case [68]. This result has been used to check our 
code for the simulation. 
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Figure 1: Various observables are plotted against a = a^/N for N = 4,8,16. The dashed lines 
represent the one-loop results at large N. 

By performing one-loop calculation around the fuzzy sphere Ajj, , we obtain the lead- 
ing large- iV behaviors as (See appendix C for the details) 

ltr(^) 2 ^ia 4 + 6, (3.4) 

1 1 Q 

(AT) ~ --a 3 + , (3.5) 



N 6 

^(>(V)^ 2 + 0, (3-6) 
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where we have introduced the rescaled parameter 

a = ay/N . (3.7) 

In the r.h.s. of eqs. (3.4) ~ (3.6) the first term represents the classical result, and the second 
term represents the one-loop correction. 




Figure 2: The history of (jj tr (A fl ) 2 )/N is shown for various a at N = 4. The horizontal lines 
represent the one- loop results for (^tr (A^) 2 ). 

In figure 1 we plot the results for (-^ tr (-F^) 2 ) and (M) obtained by Monte Carlo 
simulations. We find that Monte Carlo data agree with the one-loop results even at a = 0. 
This is rather surprising since the expansion parameter in the perturbative calculation (at 
finite N) is As we will see shortly, however, the system actually changes its behavior 
at a cc ^t=, and the agreement in figure 1 below that point should rather be considered as 
accidental. 

Let us then consider the quantity (^tr(A M ) 2 ), which we have postponed since it 
involves a subtle issue. At a = this quantity is actually divergent even for finite N, 
as first observed in numerical studies [79] and further confirmed by ref. [14]. (Ref. [77] 
provides some analytical explanation.) On the other hand, from perturbative calculations 
around the fuzzy sphere, we obtain a finite result for finite iV (See eqs. (C.6) and (C.7)). 
In order to clarify the situation, let us first look at the history of -^tr( J 4 /i ) 2 , which is 
plotted in figure 2 for various a at iV = 4. The horizontal axis represents the number of 
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"trajectories" in the hybrid Monte Carlo algorithm (See appendix A.) with the parameters 
v = 100 and At = 0.01. For a = the history has a lot of spikes, and these spikes are 
responsible for the divergence of (jjtr(A^) 2 ). As we increase a the spikes become less 
frequent, and their height gets lowered. At a > 1.1 the history looks quite regular. We 
also looked at the history at larger N, and find that it becomes regular (no spikes) for 
a > 0.5 at N = 8 and for a > 0.3 at N = 16. The transition point a tr is roughly consistent 
with a t r oc jj (i.e., a tl oc -^=). 

According to ref. [78], switching on a does not change the convergence properties of 
the matrix integrals. Therefore, unless there is some special mechanism for canceling the 
leading divergence, we get tr {A^) 2 ) = oo even for a^0. The reason why we get finite 
results by perturbative calculations should then be that such a calculation only include the 
region of the configuration space near the fuzzy sphere solution so that it does not take 
into account the configurations that have large tr (yl^) 2 , which may actually contribute 
crucially to the vacuum expectation value of that quantity. 




Figure 3: Monte Carlo results for ^= are plotted against a (left) for TV = 4, 8, 16. On the right 
the same data are plotted in a different scale. The straight lines represent the classical result 
R = | a VN for the fuzzy sphere solution. 

One can, however, define a finite quantity [14] 

R=(^±tr(AJ*y (3.8) 

which is finite at a = and behaves as 0(1) at large N [11] 3 in the present parametrization 
of the action. Due to the argument of ref. [78], this quantity is finite also for a^0. We 
plot the Monte Carlo results for R/y/N in figure 3 on the left as a function of a for N = 4, 
8 and 16. At large a the data agree very well with the classical result R = \ a y/N for the 
fuzzy sphere solution. 

In figure 3 on the right we plot R against a. We expect that R approaches a finite 
value in the large-A limit for each a in the small-a regime, and our data are roughly 

3 Strictly speaking, ref. [11] studies a slightly different observable, but the large- iV behavior is expected 
to be qualitatively the same. 
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consistent with this picture. If we assume the transition to take place at the point where 
the fuzzy sphere result R ~ ^ a y/N becomes comparable to the pure super Yang-Mills 
behavior R ~ 0(1), the transition point should be a oc which roughly agrees with the 
point a tr , where the spikes in figure 2 get suppressed. 

Thus we conclude that the system actually undergoes some transition at a oc 
below which the system behaves similarly to the pure super Yang-Mills model (a = 0). 
The quantities in figure 1 are insensitive to the transition since their behavior in the pure 
super Yang-Mills phase just happens to be the same as in the fuzzy sphere phase. Note that 
(M) a= Q = due to parity symmetry Ai ^ — -Aj, whereas from the perturbative expansion 
around the fuzzy sphere, one finds that the one-loop contribution to (M) is absent due 
to supersymmetry 4 , and as a consequence the result (C.9) has a smooth extrapolation to 
a = 0, which agrees with (M) a= o = 0. Since (M) and (^trF 2 ) are related to each other 
through the exact result (3.3), the agreement of the former propagates to that of the latter. 
Thus the success of one-loop results for these quantities in the small-a regime does not 
mean that we are still in the fuzzy sphere phase, but it simply means that these quantities 
are insensitive to the transition from the fuzzy sphere phase to the pure super Yang-Mills 
phase. 

The above results are in striking contrast to those obtained in the bosonic model [68], 
where we observed that the fuzzy sphere becomes unstable at some finite a and various 
quantities show a hysteresis behavior, implying a first order phase transition. 

3.2 Theoretical understanding based on the effective action 

In the previous section we observed that the fuzzy sphere is stable in the large-iV limit 
at any finite a unlike the bosonic model. Here we would like to provide some theoretical 
understanding of the striking difference between the bosonic and supersymmetric cases 
based on the one-loop effective action. For that purpose let us consider a one-parameter 
family of configurations given by 

A i = f3L { f ) fori = 1,2,3, 

^4 = 0, [ ' ) 

where the fuzzy sphere solution (2.6) corresponds to (3 = a. The one-loop effective action 
around (3.9) can be calculated along the line described in appendix B, and we get the 
result at large TV as 

^ <lo P = Q /3 4 " \& /3 3 ) " log N , (3.10) 

where (3 = (3 y/N. The one- loop effective action has a minimum at (3 = a for arbitrary a. 

In analogous calculations in the bosonic models [68,70,73], the one-loop contribution 
gives rise to a term proportional to log/3. Due to this term the (local) minimum disap- 
pears below some critical a, which indeed agrees well with the Monte Carlo results. In the 

4 Note that the cubic term, which breaks supersymmetry softly, does not appear in the relevant one-loop 
calculation. This is not the case, however, at higher loop calculations. 
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present supersymmetric case, the /3-dependent one-loop term is absent due to supersym- 
metry. Thus we can understand the qualitative difference between the bosonic case and 
the supersymmetric case observed in Monte Carlo simulations. 

4. Geometrical structure 

In this section we study the geometrical structure of the dominant configurations in the 
supersymmetric model. For that purpose we consider the "Casimir operator" 

Q = (4-1) 

and define its eigenvalue distribution f(x) as 



1 N 



j'=i 



where Xj (j = 1,2, ••• , N) represent the eigenvalues of Q. Let us note that (77 tr {A^) 2 ) 
discussed in the previous section is related to f(x) as 



(^tr(^) 2 ^ = (^tr(^ = {Lj^ A ^ = j\f(x)dx . (4.3) 

In figure 4 we plot the eigenvalue distribution for the same set of a and N as in figure 2. 
4.1 Power-law tail 

The results for a = reproduce the power-law behavior 

f(x) cc x- 2 , (4.4) 

at large x, which has been first discovered in ref. [79] and studied also in ref. [80]. This is 
related to the divergence of (-^ tr (A^) 2 ) discussed in section 3.1. In fact this quantity is 
known to diverge logarithmically [79], which explains the power in (4.4) due to (4.3). 

At a = 0.7 the magnitude of the power-law tail becomes much weaker, but our data 
are consistent with the existence of the power-law tail. At a > 1.1 the power-law tail 
becomes hardly visible, which corresponds to the disappearance of the spikes in the history 
of jf tr (A^) 2 seen in Fig. 2. If we assume that tr (^) 2 ) = 00 for a / 0, as we argued in 
the previous section, the power-law tail should be still there, but simply hidden by the main 
contribution coming from the fuzzy-sphere- like configurations 5 . Considering that those 
configurations are enhanced by the Boltzmann weight e const - a4N2 a t l ar g e a compared with 
the configurations that give the power-law tail, we expect that the magnitude of the power- 
law tail decreases as e - const - Q4Ar2 . Therefore the power-law tail is expected to disappear 
completely if we take the large- N limit at fixed a. 



5 Within perturbative expansion around the fuzzy sphere configuration, we expect to obtain a distribution 
which decays faster since the n-th moment J dxx n f(x) = j^(tr {(A^) 2 }™} can be calculated as a finite 
quantity to all orders. 
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Figure 4: The plot of the eigenvalue distribution f(x) of the Casimir operator Q for various a at 
N = 4. The dashed lines in the two upper figures represent a fit to the power-law behavior (4.4), 
and the vertical lines in the two lower figures represent the classical results (^-function) for the 
single fuzzy sphere solution. 

4.2 Spherical geometry 

In order to clarify the geometrical structure of the fuzzy-sphere- like configurations, we 
decompose the Casimir operator Q as Q = Q( 123 ) + Q^ 4 \ where 

Q (123) = E (A t f , (4.5) 
i=i 

Q {4) = (A,) 2 , (4.6) 

and calculate the eigenvalue distribution for Q^ 123 ) and which we denote as /^ 123 ^(x) 
and /( 4 )(x), respectively. Figure 5 shows the result for a = 1.5 at N = 4. 

The figure on the right shows that the eigenvalues of is quite small, which can 
be understood from the one-loop result (C.7) for (-^ tr (^4) 2 ) 1 _ loop ) which vanishes as 
0(-^>z log N) in the large- N limit with fixed a. As a consequence the distribution /( 123 )(x) 
shown on the left of the figure 5 is almost identical to f(x) shown on the bottom right of 
figure 4. We also find that /( 123 ^(x) is peaked around the classical result. Thus we confirm 
that the dominant configurations indeed have the geometry of a 2-sphere. 
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Figure 5: The functions /( 123 ) (x) and /( 4 ) (a;), which are the eigenvalue distributions of the operator 
q( 12 3) an( j Q( 4 ) 5 respectively, are plotted for a = 1.5 at N = 4. The classical result for the single 
fuzzy sphere solution is represented by the vertical line (<5-function) in the left figure. 

At a = the distribution f{x) has an empty region around x = 0. Similar behavior has 
been observed also in the bosonic model [68], and it can be understood by the uncertainty 
principle. Therefore the geometrical structure of the dominant configurations at a = 
should rather be considered as that of a solid ball. 



5. Dynamical gauge group 

In fact the equation of motion (2.5) has a class of solutions of the form 

Af S2) = a {4 n) l fc ) fori = 1,2,3, 

4 fcs2 > = o , 

where N = nk. These solutions represent the k coincident fuzzy S 2 , and the Casimir 
operator Q takes the value 

Q = - A {n 2 -l)a 2 l N , (5.2) 

meaning that the "radius" of the fuzzy spheres is given by p = \fn 2 — 1 ~ ^ a N, which 
becomes smaller as k increases. By expanding the theory around such a configuration, one 
obtains noncommutative Yang-Mills theory with the U(fc) gauge group [23]. 

In order to study the stability of such a configuration, we perform Monte Carlo sim- 
ulation for N = 8 and a = 1.0 using the k = 2 coincident fuzzy spheres as the initial 
configuration. In figure 6 we plot the history of the eigenvalues of Q. The horizontal axis 
represents the number of "trajectories" in the hybrid Monte Carlo algorithm (See appendix 
A.) with the parameters v = 100 and At = 0.0001. Thus the k = 2 multi fuzzy sphere 
configuration is unstable and decays into the single fuzzy sphere. 



(5.1) 



This phenomenon can be understood 
by considering the free energy, which is 
calculated in appendix B up to one-loop. 
At large N the result reads 

— W$ k ] = — a 4 -logiV 

N 2 !- 1oo p 24 k 2 & 

Since the one-loop term represented by 
the second term is independent of k, we 
find that the free energy takes the small- 
est value for k = 1. 

Although the conclusion concerning 
the dynamical gauge group is the same as 
in the bosonic models [68,70,73], we note 
that the reasoning is different. In the 
bosonic models the one-loop term in the 
free energy has the — log k 2 term, which 
actually favors large k. However, if one 
decreases a so that the one- loop term becomes more important, the fuzzy sphere solutions 
disappear before the free energy starts to favor k > 1. In the present supersymmetric case, 
the fuzzy sphere solutions remain to be there, but due to the absence of the — log k 2 term, 
the k = 1 solution is always favored. 

6. Conclusion 

In this paper we have studied the dimensionally reduced 4d super Yang-Mills model with 
an extra Chern-Simons term, which incorporates fuzzy spheres as classical solutions. We 
have found that supersymmetry indeed has substantial effects on the dynamics of fuzzy 
spheres. 

While the observables that appear in the action change continuously as we vary a, the 
model actually possesses two distinctive phases, which is demonstrated by a well-defined 
observable ^ \J tr (A^) 2 ^ . The tail of the eigenvalue distribution changes drastically as 
one moves from one phase to the other. In the pure super Yang-Mills phase we observe the 
same power-law tail as the one known for a = 0, but it disappears in the fuzzy sphere phase 
in the large- N limit. This allows us to identify the critical point quite accurately. From 
our Monte Carlo data up to N = 16 and some theoretical considerations, we speculate that 
the transition point a tr goes to zero as <Str = O • 

Our results are in sharp contrast to the results obtained for analogous bosonic models 
[68, 70, 73], where the fuzzy sphere becomes unstable at some finite critical a. A strong 
first-order phase transition has been observed in the bosonic models, and the (lower) critical 
point obtained by Monte Carlo simulation can be reproduced very well from the one-loop 
effective action. In the present paper we have shown that the one-loop term in the effective 
action vanishes up to an irrelevant constant due to supersymmetry. This explains our 
observation that the fuzzy sphere is stable down to vanishingly small a at large N. 
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Figure 6: The history of the eigenvalues of the 
Casimir operator Q obtained for N = 8, a = 1.0 
using the k = 2 solution as the initial configura- 
tion. The horizontal lines represent classical results 
for the single fuzzy sphere (k = 1) and the two co- 
incident fuzzy-spheres (k = 2), respectively. 
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We have also studied the dynamical generation of the gauge group. In the bosonic 
case the quantum instability of the fuzzy spheres was an obstacle in obtaining a non-trivial 
gauge group in the true vacuum. In the present super symmetric case, this instability is 
gone, but we also lost the one-loop term in the free energy for the coincident fuzzy spheres, 
which favors higher multiplicity. As a result, we obtain the U(l) gauge group again. We 
note, however, that this conclusion applies only to the fuzzy sphere phase, and whether 
we can obtain a nontrivial gauge group in super symmetric models without a Chern-Simons 
term such as the IIB matrix model [2] is still an interesting open question. 
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A. Details of the Monte Carlo simulation 

In this section we explain the algorithm used for our simulation. Our algorithm is similar 
to the one adopted in ref. [11], but the crucial difference is that we make the Metropolis 
reject /accept procedure at the end of each trajectory. In the previous algorithm there was 
a systematic error due to discretization required for solving the Hamiltonian equation, and 
the step size At for the discretization had to be sent to zero. In the present algorithm 
we do not need such an extrapolation. Another difference is that we do not use the noisy 
estimator for estimating the r.h.s. of the Hamilton equation since it causes some systematic 
error. Instead we invert the Dirac operator directly using the LU decomposition. Each 
of the two modifications increases the computational effort for making one trajectory (for 
fixed parameters in the algorithm) from 0(iV 5 ) to 0(N 6 ), which is the price we have to 
pay to make the algorithm "exact" . 

The "exact" algorithm is essentially the hybrid Monte Carlo (HMC) algorithm, which 
is used in studying the large- N behavior of the phase quenched version of the IIB ma- 
trix model [12]. In that case the one-loop approximation has been used to decrease the 
computational effort from 0(iV 6 ) to 0(iV 3 ). The hybrid algorithms are standard in full 
QCD simulations, and it is useful also in simulating matrix models as demonstrated in 
refs. [11,12]. If we had used the Metropolis algorithm [14] in the present model, for in- 
stance, the computational effort would have been 0(N 8 ). Note, however, that simulating 
matrix models is generally harder than simulating field theories due to the non-local nature 
of the interaction. Even in the bosonic case, the computational effort is at least 0(N S ), 
which grows faster than the number of d.o.f., which is 0(N 2 ). 
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Let us first recall an explicit form of the fermion determinant derived in ref. [11]. We 
define a complete basis for the general complex N x N matrices as 



(t a )ij = SaJ jja (a = 1,2,-.. ,N 2 ) , (A.l) 

where i a ,ja are integers within the range 1 < i a ,ja < N satisfying 

a = N(i a -l)+j a . (A.2) 

By taking into account that the fermionic matrices ip a and ip a are traceless, integration 
over the fermionic matrices yields the fermion determinant det M., where the 2{N 2 — 1) x 
2{N 2 — 1) matrix Ai is given by 

M aaW = M' actjbp - M' N 2 a bf3 S iaja - M' aa N 2p S ibjb . (A.3) 

Here the 2N 2 x 2N 2 matrix M' is defined as 

M'^^iT.Utr^lA^}). (A.4) 

The effective action for can be written as 

S cS [A] = S h [A] - log Det M [A] . (A.5) 

Following the idea of the hybrid Monte Carlo algorithm, we introduce auxiliary bosonic 
hermitian matrices P M and consider the action 

5hmc[P, A] = l - tr (i^) 2 + S cS [A] . (A.6) 

Since P p does not couple to A^, we retrieve the original model trivially by integrating out 
P M . We regard the action SmAc[P,A] as the Hamiltonian of a classical system described 
by A^{t) and its conjugate momentum P M (r), where r denotes the fictitious time of the 
classical system. Then as an update procedure, we may take the old configuration (P, A) 
as the initial configuration (-P(O), A(0)) and solve the Hamiltonian equation 

dr - 0(P„) y - ™< ' (A - 7) 

d(Pfi)ij _ dSftMC _ dS e ff 



dr diA^ij d{A ll ) ij 

= n(-[A v ,[A^A v }} + 2iae^ vp A u A p ) - Tr (AT 1 -^- 



(A.8) 



for a finite fictitious time T (this defines "one trajectory") to obtain (P(T), A(T)). The 
symbol Tr in (A.8) denotes a trace over the 2 (N 2 — l)-dimensional index, and the derivative 
g(X^ is given explicitly by 
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Since the trace of is not conserved during the evolution, we subtract the trace part 
A'^ = A^{T) — {-^ tr ^4 M (T)} 1. Thus we obtain the updated configuration (P',A r ), where 
P'^ = Pfj,(T). Due to the Hamiltonian conservation, this update procedure preserves the 
action Sjimc[P, A]. Using also the fact that the transition between (P, A) and (— P', A') 
is reversible, one can readily verify the detailed balance. After each trajectory, we update 
the momentum P^ fixing A„, which can be done by simply generating gaussian variables 
since the P^-dependent part of the action (A. 6) is gaussian. This procedure is necessary 
to avoid the ergodicity problem. 

In actual calculations we have to discretize the Hamiltonian equation (A. 8). The re- 
versibility of the time evolution can be preserved by using the so-called leap-frog discretiza- 
tion, but the Hamiltonian conservation is inevitably violated. However, we may accept the 
configuration (P',A') as the updated configuration with the probability max(l, e~ ASliMC ), 
where AShmc = "Shmc[-F"> A'\ — Shmc[P ; ^4]) an d duplicate the old configuration when 
rejected. By adding such a Metropolis accept/reject procedure, we can preserve the de- 
tailed balance. The step size At for the time evolution should be small enough to keep the 
acceptance rate reasonably high. The discretized Hamiltonian equation is given by 

(4 1) )^ = (4 0) )^ + Ar(P / l 1 /2) ) .. ) (A . 10) 

{p { r 1/2) h = (pt l/2) h - A -^-(4 n) ) > 

(4 n+1) )^- = + At ( p !? +1/2) h > (A-ii) 
(p^h = (Pt 1/2) h - ^ > ( A - 12 ) 

where n = 1,2, 1 and T = u At, and we have introduced the short-hand notation 

pff^ = P M (rAr) and A^ = A M (r At). At each step of the "Molecular Dynamics", we 
have to calculate the inverse M -1 , and at the end of each trajectory, we have to calculate 
detM. These are the dominant part of the numerical calculation, and it requires a CPU 
time of the order of 0(A^ 6 ). 

The hybrid Monte Carlo algorithm involves two parameters T and At, which can 
be optimized in such a way that the computational effort for obtaining one statistically 
independent configuration is minimized. The optimization can be done in a standard 
way [12]. First we fix T and optimize At so that the effective speed of motion in the 
configuration space, which is given by the acceptance rate times At, is maximized. Using 
the At optimized for each T, we minimize the autocorrelation time (in units of "Molecular 
Dynamics step") with respect to T. For instance, at iV = 16 and a = 0.0 we obtain the 
optimal values At ~ 0.006 and T ~ 1.0. 

B. One-loop free energy 

In this section we formulate the perturbation theory around fuzzy sphere solutions, and 
derive the one-loop free energy. We decompose A^, -ip and V> into the classical background 
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and the fluctuation as 

A M = X M + i M , (B.l) 
V> = x + ^, ?P = x + 4>, (B.2) 

and obtain the free energy around the classical solutions by integrating over A^, ip and 
V> perturbatively. Here we take the classical solution to be the k coincident fuzzy spheres 
Xn = \ x = X = 0, which includes the single fuzzy sphere as a special case k = 1. 

In order to remove the zero modes associated with the SU(iV) invariance, we introduce 
the gauge fixing term and the corresponding ghost term 

S gS =-^tr[X, 1 A,f, (B.3) 
S ghost = -Ntr{[X^c][A^c]) = Nte{c[X^[A^c]]) , (B.4) 

where c and c are the ghost and anti-ghost fields respectively. 

The total action S to tai = S + S g j, + S g \ lost can be written as 

Stotal = <Scl + "Skin + Shit , (B-5) 

S cl = Ntr I ~ \X^X v f + f *<* £ ^ Xi Xj X k J , (B.6) 
S kin = N tr - [Ap, A v ] X v ] + i a ^ e ijk [ A h Aj } X k 

\ i,j,k=l 
- l -[X^A u ] 2 + c[X^ [X„, cW-^T^X^ , (B.7) 

S- m t = Ntr {-[A^A V ][X^A U ] - i [A^A V \[A^A V } 

2 3 \ 

+ 3«a CijkAiAjAk + ciX^, c]] -^r M [i M , 4>] \ . (B.8) 

jj,fe=i / 

The linear terms in A^ cancel since X^ satisfies the classical equation of motion. 

Noting that the background configuration X^ includes a factor of a, we can rescale 
the fluctuations as An ^ ctA^, c i— > a c, c i— > a c, V l— ^ a5 V'j V' l— ^ a ^ V' so that all the 
terms in the total action Stotal become proportional to a 4 . This means that the expansion 
parameter of the present perturbation theory is ^ . 

The free energy W is defined by 



e -w 



= J didcdcd^d^e" 5 ^ 1 , (B.9) 

which can be calculated as a perturbative expansion W = Yl'jLo ^0 > wnere Wj = 
0(a 4< - 1- ^) . The classical part is simply given by Wo = S^[X] . In order to evaluate the 
one- loop term W\, we note that the kinetic terms can be written as 

Skin = Ntr (^A u (V x ) 2 A u + c{V x ) 2 ^ -Ntr (^T^P^) , (B.10) 
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where we have introduced the operator P M which acts on a traceless N x N matrix M as 

^M=[X M ,M]. (B.ll) 

Then the one-loop term can be expressed as 

W 1 = W 1>b + W lf , (B.12) 
W 1>h = Tr log {N{V^) 2 } , (B.13) 
W hi = -Tr' log (NT^) , (B.14) 

where the symbol Tr denotes the trace in the (iV 2 — l)-dimensional linear space which 
consists of traceless N x N matrices, and Tr' includes the trace over spinor indices as well. 

B.l Single fuzzy sphere 

Let us first consider the single fuzzy sphere = . The classical part is given by 

W = -j l N 2 a 4 (N 2 -l) , (B.15) 

and the one-loop terms can be written as 

W 1>h = Tr log (TV a 2 Q) , (B.16) 
W hi = -Tr' log ( N a V ) . (B. 17) 

The operators Q and V are defined as 

3 3 

Q = Y J (d) 2 , V = Y J «i£i, (B.18) 

i=l i=l 

where L% acts on a traceless N x N matrix M as C % M = [L^ , M] . In order to evaluate 
the one-loop terms, we need to solve the eigenvalue problem of the operators Q and V. 

The eigenvectors of the operator Q are given by the "matrix spherical harmonics" Yj m 
(I = 0, 1, ■ ■ ■ ,N—1 and m = —I, - ■ ■ ,1), which span a complete basis of the space of N x N 
matrices and have the properties analogous to the usual spherical harmonics such as 

^ tr i^lmXl'm') = $ll'8mm' , (B.19) 

YL = (-l) m ^,- m • (B.20) 
The corresponding eigenvalues are given by I (I + 1); i.e., 

QY lm = l(l + l)Y lm . (B.21) 

Thus the one-loop term from the bosonic contribution is obtained as 

JV-l 

Wi, b = (2Z + l)log[iVa 2 Z(Z + l)] . (B.22) 
l=i 
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Here I = has been omitted from the sum since the trace Tr in (B.13) should be taken in 
the space of traceless N x N matrices. 

In order to solve the eigenvalue problem of the operator V, we note that 

V = ^m 2 -Q~l, (B.23) 
i=i 

where we have defined the "total angular momentum" operator 

Ji = A + f ■ (B.24) 

By making a linear combination of eigenvectors of Q with the eigenvalue 1(1 + 1), we can 
construct the eigenvectors of both Yli=i (-^) 2 an d Jz with the eigenvalues j (j + 1) and m, 
respectively, where j can be either j = 1+^ (Z = 0, ■ ■ ■ , N— 1) or j = I— \ (Z = 1, ■ • ■ , iV— 1), 
and m takes half-integer values in the range \m\ < j. Explicitly, the eigenvectors are given 
by the "matrix spinorial-spherical harmonics" 



I + - + m I + - — m 



y\-l, m = y ^TTi-^ ® IT) - y ^f^^, m+ ^ ® 1 1 ) . (B-26) 

for the cases j = Z + ^ and j = I — ^, respectively. Here the symbol 1 1 ) and | { ) denotes the 
two-dimensional eigenvectors of 03 corresponding to the eigenvalues 1 and —1, respectively. 
From eq. (B.23), the "matrix spinorial-spherical harmonics" are also eigenvectors of V with 
the eigenvalues 

V = j{j + l)-l{l + l)-- A (B.27) 



I for j = I + i 

-{1 + 1) for j = /-i. 



(B.28) 



Namely we have the relation 

vy l+hm = iy l+hm , (B.29) 
vy' l _ lm = -(i + i)y' l _ hm . (B.30) 

Thus the one-loop term from the fermionic contribution is obtained as 

nv-i jv-i 

2(Z + 1) log (iVaZ)+ 21 log{iVa(Z + l)} 



W 



i,f 



1=1 1=1 



(B.31) 



where Z = has been omitted from the first sum since the trace Tr' in (B.14) should be 
taken in the space of traceless N x N matrices (and over the spinor indices) . Let us rewrite 
the above expression as 

W hi = -W hh - {N 2 - 1) log N + log N . (B.32) 
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Thus the fermionic contribution cancel the bosonic contribution up to the a-independent 
constant. From (B.15) and (B.32) the one-loop free energy for the single fuzzy sphere is 
obtained at large N as 

w^ loop ~ n 2 (__L „4 _ logN ^ _ (B 33) 

B.2 k coincident fuzzy spheres 

Next we consider the k coincident fuzzy spheres = . The classical part of the free 

energy is given by 

W = -^a i (n 2 -l) . (B.34) 

In order to calculate the one-loop term, we consider the n x n version of the matrix 
spherical harmonics Y^ , and define 

Yti b) - ® * {a ' b) , (B-35) 

where e^ a ' 6 ^ denotes a k x k matrix whose (a, b) element is 1 and all the other elements are 
zero. The N x N matrices Y^f^ form a complete basis of the space of N x N matrices, 
and they are the eigenvectors of the operator (V^) 2 for the present background; i.e., 

{^fY^=a 2 l{l + l)Y^ . (B.36) 

Let us note that the operator (P^) 2 has k 2 zero modes corresponding to / = m = with 
arbitrary (a,b). The mode Yla=i ^oo'°^ corresponds to the trace mode, which should be 
omitted due to the traceless condition. Here we omit the other zero modes by hand, 
and simply consider the non-zero modes. Then the one-loop term Wi^ from the bosonic 
contribution is obtained as 

n—l 

^i,b = fe 2 ^(2/ + l) log[iVa 2 /G + l)] . (B.37) 
1=1 

The calculation of the fermionic one-loop term proceeds in the same way except that 
we have to use the matrix spinorial-spherical harmonics for each of k 2 blocks. In this case 
we have 2k 2 zero modes, and two of them correspond to the trace mode. We omit the 
other zero modes by hand and obtain 



n— 1 n—l 



2{l + l)log(Nal) + Y J 2/ log{ Na{l + 1)} 



i=i 



(B.38) 



The cancellation between the bosonic contribution and the fermionic one occurs here as 
well, and the one-loop free energy is given by 

^l-loop - iV 2 ("2^2 « 4 - l0 § N ) ■ ( B - 39 ) 
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C. One-loop calculation of various observables 

In this section we apply the perturbation theory discussed in the previous section to the 
one-loop calculation of various observables which are studied by Monte Carlo simulations 
in this paper. We take the background to be k coincident fuzzy spheres X^ = , but 

the results for the single fuzzy sphere can be readily obtained by setting k = 1. As in 
appendix B.2, we omit the zero modes for k > 2. 

We note that the number of loops in the relevant diagrams can be less than the order of 
1/a 4 in the perturbative expansion since we are expanding the theory around a nontrivial 
background. At the one-loop level, the only nontrivial task is to evaluate the tadpole 
((A^)ij) explicitly, which, however, turns out to vanish due to supersymmetry. 

C.l Propagators and the tadpole 

The propagators for A^, the ghosts and the fermion fields are given respectively as 

ab 1=1 m=-l \ ' / J 

_. n— 1 I i 

(<cW(a>„) = i££ £ w^TTT) ( y -% ( y -"% ' < C ' 2 > 

ab 1=1 m=—l 

^EE m |^(^„)>^) w .(c,) 

where the symbol ( • )o refers to the expectation value calculated using the kinetic term 
Skm only. The tadpole (^4j)i-i 00 p (i = 1>2,3) at the one-loop level can be calculated as 

(A)i-ioop = (NAi tr ([A v ,A p ][X,,,A p ])} q - (nA^v (c[X v , MUc]])) q 

— (jSlAi tr ($T v [A v ,^) q . (C.4) 

By redoing the calculation in ref. [35] in the present model, we find that the bosonic 
contribution and the fermionic contribution cancel each other even at finite N . We also 
find that (^4} = to all orders in perturbation theory due to parity invariance A4 1— > —A4. 

C.2 One-loop results for various observables 

Using the propagator and the tadpole obtained in the previous section, we can evalu- 
ate various observables easily at the one-loop level. For instance the two-point function 
(^tr(A^) 2 ) can be evaluated as follows. Let us decompose it as 

(k tr < - 4 " )2 } = (71 £ tr (A>2 ) + (k b {A » 2 ) ■ (&5) 
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Each term on the r.h.s. can be calculated as 



N 



i=i 



1— loop 



1 3 

- [tr(XiXi) + 2tr (x i {A t ) 1 _ loop ) + (tr(i 4 ) 2 ) 



i=l 

1/2 \ 3 ^ 2/ + 1 

4 ( - - 1) + + ^g7(7TT) 



]v tr(i4)2 



1— loop 

At large iV with fixed a = a y/~N, we get 



Q" 



1(1+1) 



(C.6) 
(C.7) 



1 / 1 



N \N 



tr (^) 



11 

N \ N 



i=i 



1— loop 

The Chern-Simons term (M) can be evaluated as 

2z 



1— loop 



(C.8) 



z % r / ~ \ 

(M)i_i oop = ~^ e ijk tr{XiXjX k )+ZtV \ XiXj{A k )i-\oopj 



--a*(n 2 -l). 



At large iV with fixed a = «yiV, we obtain 

1 



(-^Ol-loop ~-6jfe2«' 



(C.9) 



(CIO) 



The observable (j^trF 2 ) can be calculated in a similar manner, but we can also obtain it 
from the exact result (3.3) using (C.10) as 



T7tr(i^) 2 



1— loop 



6^1-^ -3a (iWVioop (C.ll) 

(C.12) 



~ Q; + 6 



C.3 Alternative derivation 

Since trF 2 and M are the operators that appear in the action S, we can obtain their 
expectation values easily by using the free energy calculated for the k coincident fuzzy 
sphere in Appendix B. Let us deform the bosonic action as 



S h {/h,P2,a) = N 2 



1 



0! ti(F^) 2 + (3 2 aM 



with two free parameters fii, 02, and define the corresponding free energy as 



(C.13) 



(C.14) 
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Then (tr(F^) 2 ) and (M) can be obtained by 

(M) 

_1 1 1 

By rescaling the integration variables as i-> j3 1 4 A^, if; (3% ij) and ij) i-> /3f Y>, we get 

W(/3i,/? 2 ,«) = ^(iV 2 -l) log ft + W(l,l,a/3r ! /3 2 ) • (C.17) 
Using the one-loop result 

iV 2 

W(l, 1, a)i_i oop = -— a 4 (n 2 - 1) - fe 2 (n 2 - 1) logiV + k 2 logn , (C.18) 
we can reproduce (C.9) and (C.ll). 
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